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Abstract: Salt expansion in sulfate saline soils that are widely distributed in northwestern China causes 
serious infrastructural damages under low-temperature conditions. However, the mechanism of salt 
expansion under low temperatures is not clear. In this study, we conducted a series of cooling experiments 
combined with salt crystallization to study this mechanism, and employed an ionic model to calculate the 
supersaturation ratio of the solution. During the experiments, the strength and the process of salt 
expansion were examined under different cooling rates and various crystal morphologies. The relationship 
between temperature and supersaturation ratio under transient conditions was also considered. Results 
indicate that the initial supersaturation ratio of a sodium sulfate solution is closely related to 
environmental conditions, and that this ratio decreases with slowing the cooling rates and stabilizing the 
crystal forms. Higher initial supersaturation ratios lead to an increased non-steady-state zone, resulting in 
less salt expansion. On the other hand, chloride ion content has a distinct influence on the crystallization 
supersaturation ratio of the sodium sulfate solution, and higher chloride ion content can inhibit salt 
expansion in sodium saline soils. These findings help explain salt expansion mechanisms in complex 
conditions such as seasonally frozen soils, and thus help search for improved methods of preventing salt 
expansion in sulfate saline soils. 
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1 Introduction 


Saline soils are widely distributed in north western regions of China and cover 2% of the total land 
area (Luo, 1980). Salt expansion, associated with the chemical and physical properties of saline 
soils, can seriously damage infrastructures such as roads, canals, and airports. However, with 
increasing economic development in western China, these infrastructures are set up densely in 
saline soil regions, thus greatly elevating the social and economic risks to the salt expansion. 
Therefore, study of salt expansion in sulfate saline soil is needed for reducing those social and 
economic risks. 

Previous studies of salt expansion in sodium sulfate soils in China mainly focused on 
measurement of the expansion rates of saline soils under different salt contents, water contents, 
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and other environmental conditions (Li et al., 1989; Yuan and Li, 1995; Xu et al., 1995, 2010; Wu 
et al., 2001). Recently, predictive models for salt expansion in sodium saline soils were developed 
from experimental studies (Gao et al., 1996; Chen et al., 2006; Zhao, 2012). And the calculations 
of salt expansion and precipitation of salt crystals in soils were usually depended on the phase 
diagrams of the solution in a steady-state. However, these predictive models are unable to 
accurately estimate salt expansion in a complex system under natural conditions. Generally, the 
fundamental cause of salt expansion is salt crystallization, and supersaturation is a 
thermodynamic status for salt crystallization. That is, when the supersaturation ratio of the salt 
solution exceeds initial supersaturation ratio, salt crystallization begins and salt heaving of the 
soil occurs (Mullin, 1993; Mersmann, 2001; Zhang, 2011). The Pitzer model has been widely 
used to calculate supersaturation in porous media (Steiger, 2005; Espinosa et al., 2008), but most 
related studies just considered salt crystallization and the initial supersaturation ratio in brick, 
concrete, rock and other construction materials, instead of soils (Mersmann, 1995, 2001; Rocha 
and Martina, 2002; Espinosa et al., 2008; Derluyn, 2012). When soil was examined, salt 
expansion was only studied at room temperature (Mokni et al., 2010). Therefore, salt expansion at 
lower temperatures and under non-steady-state conditions remains uninvestigated, and damage to 
infrastructure caused by salt expansion remains troublesome in cold regions. 

Researches on sulfate saline soils in cold regions must consider the following questions: At 
what supersaturation ratio does salt crystallization occur in soils at low temperature? How does 
the initial supersaturation ratio of solutions in soil vary under different environmental conditions? 
How does salt crystallization influence salt expansion? Therefore, given the large variations in 
temperature on the Tibetan Plateau, this study aims: (1) to explore the crystallization status of 
sulfate under non steady-state conditions; (2) to identify the relationship between crystal 
morphology and salt expansion; and (3) to analyze the influence of chloride ion content on the 
supersaturation ratio of sulfate saline soils. The results may be useful for understanding the 
properties of saline soil, and in studying land surface processes in cold regions with sulfate saline 
soils. 


2 Materials and methods 


2.1 Calculation of supersaturation ratio 


The fundamental cause of salt expansion is crystallization, as salt heaving cannot occur without 
the presence of salt crystals. However, salt crystallization is itself dependent on supersaturation as 
the necessary thermodynamic condition. Supersaturation is a metastable state during which 
crystal formation conditions shift from unsteady to steady. Crystallization occurs when the 
solution is supersaturated with respect to the crystallizing phase. If a salt contains hydrated water, 
the water-salt interaction must be considered, and the supersaturation ratio of the solution can be 
expressed as: 


U, - (£y x", (1) 


where Ua is the supersaturation ratio of salt solution; aw is the water activity; a; is the water 


activity at saturation; v is the number of ions produced per solute molecule; vo is the number of 
hydrated water molecules per salt molecule; a+ is the mean ion activity; and gf is the mean ion 


activity at saturation. The formulae for calculating these parameters are given in Appendix A. 

In saline soils, ions are mainly present in moisture, thus, the concentration of a salt in solution 
can be determined by the salt content and water content. Based on the saturation curve of a 
sodium sulfate solution (Xu et al., 1995), there will be little salt crystallization below the initial 
freezing temperature of soil, compared to that at room temperature. Therefore, we consider that 
solutes are excluded from ice crystals below freezing temperature, and the salt concentration in 
the remaining brine will increase as water freezes. The unfrozen water content is mainly 
determined by the ice-water phase change, and can be obtained by fitting experimental data (Wan 
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et al., 2017). 
2.2 Experimental design 


Sulfate soil is a predominant and representative type of sulfate saline soil found in western China 
(Wan et al., 2015), and sodium sulfate is most prevalent in large areas of the Beilu River on the 
Qinghai-Tibet Plateau. Therefore, sodium sulfate soil was selected to study salt expansion 
mechanisms. The samples consist of silty clay soils collected approximately 1 km from the 
western side of the Beilu River highway overbridge (34°54.217'N, 92?56.362'E), near the 
abandoned Qinghai-Tibet Highway. Factors that affect salt expansion include salt content, water 
content, and dry density (Xu et al., 1995; Mokni et al., 2010). Thus, water content and dry density 
used are those of the natural soil, and the influence of salt content (salt concentration) on 
expansion is studied under transient conditions. 


2.2.1 Cooling experiments using sodium sulfate soils 


Sodium sulfate soils were manufactured from the natural samples at room temperature, by adding 
the following percentages of sodium sulfate 0.5%, 0.6%, 0.7%, 0.8%, 0.9%, 1.0%, 1.1%, 1.2%, 
1.4%, 1.5%, 2.1%, 2.6%, 3.2%, and 3.8%. The water content of each soil was 18%. Thermistors 
(accurate to +0.01°C) encased in waterproof plastic sheets were inserted into the samples to 
monitor temperature variations. The test samples were put into specially-made iron boxes 
(diameter 3.3 cm, height 3.8 cm), and the dry density of the test samples was 1.84—1.86 g/cm’, 
being consistent with that of the in situ soil. The samples were placed into a cooling bath at 25°C 
and were maintained at that temperature for 30 min. Temperature of cooling bath was decreased 
from 25°C to —20 °C by a rate of 0.02°C/min, then samples were maintained at —20°C for 2 h. The 
data acquisition and data logging system was a CR3000 micrologger (Campbell Scientific 
Headquarters, Logan, UT, USA) with recording interval set to 10 seconds. The initial salt 
crystallization temperature is represented by the abrupt rise of the cooling curve (the heat released 
during salt crystallization). At negative temperatures, the initial salt crystallization temperature of 
soils with low salt content can be determined by the second phase-transition point on the cooling 
curve (Lai et al., 2016; Wan et al., 2017). 

2.2.2 Effects of cooling rate on supersaturation ratio and salt expansion 


Mirabilite is the predominant crystal type in sodium sulfate soils at high humidity, and the 
variation of the initial supersaturation ratio with various salt contents can be studied under 
different cooling rates when mirabilite crystallizes. We used sodium sulfate soils with sodium 
sulfate contents of 2.1%, 2.6%, 3.2%, and 3.8% and a water content of 18%. The samples were 
placed in plastic cups and compacted to a diameter of 6.3 cm and height of 7.5 cm. NS-WY02 
displacement sensors with an accuracy of +0.001 mm (Shanghai TM Automation Instruments Co., 
Ltd., Shanghai, China) were installed into the soil samples to monitor variations in displacement. 
A programmable ultra-low temperature chamber (Dongguan Terchy Test Equipment Co. Ltd., 
China) was used to replicate the temperature regime of the Qinghai-Tibet Plateau. Thermistors 
were inserted into the samples to monitor variations in temperature. The cooling rates were set at 
0.02, 0.1 and 1?C/min. Temperature was decreased from 25°C to 0°C, and samples were kept at 
0°C for 2 h. 

2.2.3 Effect of chloride ion content on salt expansion 

This section outlines the test employed to study the effect of chloride ion content on salt 
expansion in sodium sulfate soils. To ensure that the sulfate crystals would not crystallize at 25?C, 
we prepared three sulfate soil samples at room temperature with different salt contents. Sodium 
sulfate content was fixed at 3.8%, while sodium chloride contents were 0%, 3.8%, and 5.8%. The 
water content of these preparations was 18% (the concentration of the sodium sulfate solution was 
1.5 mol/L and the concentrations of the sodium chloride solutions were 0, 3.6 and 5.5 mol/L, 
respectively). Dry bulk density was maintained above 1.65 g/cm?. Test samples were placed into 
specially made organic glasses (diameter 6.3 cm, height 10.0 cm), and samples were then placed 
in a programmable ultra-low temperature chamber; displacement sensors were installed into the 


JOURNAL OF ARID LAND 


soil samples to monitor displacement variation. Temperature was decreased from 25°C to 0°C at 
2°C or 1°C intervals and samples were maintained at each temperature for 30 min. Soil 
displacement was measured by a date Taker CR3000. 


3 Results 


3.1 Salt expansion of sodium sulfate soil under transient conditions 


3.1.1 Initial supersaturation ratio 


The initial supersaturation ratio can be calculated using Equation 1 for mirabilite crystal 
morphology. As shown in Figure 1, the freezing point is the cut-off point for the rate change of 
the initial supersaturation ratio in sodium sulfate soils. Above the freezing point, the regression 
slope is —0.18, and the initial supersaturation ratio reaches up to 4.2. However, the initial 
supersaturation ratio increases sharply as the temperature decreases below the freezing point of 
the sodium sulfate soils. The slope of the regression line is —4.9 and the value of the initial 
supersaturation ratio reaches 45 at —11.5?C. If the supersaturated state of the solution is defined as 
an unsteady state, in which the supersaturation ratio of the solution changes from 1 to the value at 
which salt crystals begin to crystallize, then the unsteady zone of salt crystallization can be 
represented as the shaded zones in Figure 1. 
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Fig.1 Initial supersaturation ratios of solution in sodium sulfate soils 


3.1.2 Effect of cooling rate on initial supersaturation ratio 


Figure 2 shows that the initial supersaturation ratio increases as the cooling rate increases. 
Together with Figure 1, Figure 2 also shows that the area of the unsteady-state zone of the 
solution increases with increased cooling rate. In addition, the condensate depression of each 
solution, at a given temperature, increases with a faster cooling rate, thereby making initial 
crystallization increasingly difficult. It is also clear from Figure 2 that the initial supersaturation 
ratio varies significantly with cooling rate and with varying salt concentration. In the sample with 
a salt content of 3.8%, this ratio increased from 0.037 to 0.131 as the cooling rate increased 
from0.02 to 0.1°C/min and from 0.1 to 1°C/min, respectively. However, the ratio increased from 
0.226 to 0.572 for the soil with a 2.1% salt content given the same increase in cooling rates. This 
illustrates that cooling rate significantly influences the initial supersaturation ratio of the solution 
in low-salt sodium sulfate soils. 

3.1.3 Effect of cooling rate on salt expansion 


For a given salt content, the salt expansion ratio decreases as the cooling rate increases, and the 
difference in the ratio between different cooling rates increases as the initial salt content increases 
(Fig. 3). In the sample with a salt content of 2.1%, the increase in expansion ratio reaches 0.021% 
and 0.292% as the cooling rate decreases from 1 to 0.1°C/min and from 0.1 to 0.02°C/min, 
respectively. However, for sodium sulfate saline soil with a 3.8% salt content, this increase 
reaches 0.818% and 1.268%, respectively, for the same decrease in cooling rates. This also 
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demonstrates that soils with high salt content have smaller unsteady zones and that in these soils 
salt heaving is especially strengthened when the temperature drops quickly. 
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Fig.2 Initial supersaturation ratios under different cooling rates 
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Fig.3 Salt expansion ratios under different cooling rates 


3.1.4 Salt expansion process 

To highlight the effect of salt expansion, we also analyzed sodium sulfate soils with a salt content 
of 3.8%. Mirabilite is formed in cooling conditions (vo=10, v=3), and the supersaturation ratio is 
calculated by Equation 1. The changes in supersaturation ratio and salt expansion ratio with 
changing temperature are shown in Figure 4. 
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Fig. 4 Salt expansion ratios and supersaturation ratios under different cooling rates (dashed lines indicate the 
supersaturation ratio of solution and solid lines indicate the salt expansion ratio) 


When the supersaturation ratio of the sodium sulfate solution is less than 1, sodium sulfate soils 
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are in a stable state, but soil displacement decreases slightly due to soil contraction. As 
temperature decreases, the supersaturation ratio increases and exceeds 1, and the soils enter an 
unsteady state. If the supersaturation ratio is higher than the initial ratio, salts will crystallize. The 
salts crystallize first in the soil with the lowest cooling rate, leading to a salt expansion (Fig. 4). 
The increase in the cooling rate leads to an enlargement of the unsteady zone, and to an increased 
supersaturation ratio of the solution, causing the initial crystallization temperature to decrease. 
The larger unsteady zone of the salt solution suggests that the concentration of the salt solution is 
higher during the cooling process and that relatively small amounts of salt crystals are generated 
in these soils. Thus, the larger salt expansion ratios correspond to the smaller unsteady zones. 
Figure 4 also shows that the unsteady zones of salt solutions decrease continuously as temperature 
decreases following initial crystallization and that the supersaturation ratios tend towards 1. This 
also suggests that salt expansion occurs more readily when the unsteady zones of salt solutions 
contracts. 

3.1.5 Effect of sulfate sodium crystal morphology on supersaturation ratio 


The salt expansion ratio of sodium sulfate soils is directly influenced by the morphology of salt 
crystals. From Equation 1, we calculated the initial supersaturation ratios of solutions for different 
sulfate crystal morphologies (Fig. 5). Figure 5 shows that the initial supersaturation ratio 
decreases during the transformation of anhydrous sodium sulfate to sodium sulfate heptahydrate 
and sodium sulfate decahydrate. Furthermore, this reduction of the initial supersaturation ratio is 
enhanced as the temperature decreases. The initial supersaturation ratio for anhydrous sodium 
sulfate is 0.14 greater than that for sodium sulfate decahydrate at 17°C, while this difference is 
0.54 at O°C. Generally, for salt crystals to achieve stable crystal morphology, the initial 
supersaturation ratio of the solution must decrease. Thus, the process of cooling the sodium 
sulfate soils causes crystals to constantly alter from an unsteady to a steady state. 
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Fig. 5 Initial supersaturation ratio of solution with different sulfate crystal morphologies 


3.1.6 Crystal morphology under different environmental conditions 


The morphology of sodium sulfate crystals changes with environmental humidity and material 
pore size, and the growth and size of crystals are also influenced by pore size. Although mirabilite 
typically exhibits a needle-like structure (Fig. 6a), the transformation from mirabilite to 
anhydrous sodium sulfate, as the temperature reaches 32.4?C, still shows needle-like branches as 
shown in Figure 6b, whereas anhydrous sodium sulfate ultimately becomes powder-like after 
sodium sulfate soils experience wetting-drying cycles (Fig. 6c). Sodium sulfate heptahydrate has 
a pyramidal shape and is difficult to observe in soils due to its instability. The shapes of salt 
crystals in soil differ from those in solution. Typically, salt crystals in soils are very small and 
their observation requires the use of a scanning electron microscope (SEM). However, the shapes 
of some salt crystals in sodium sulfate solutions can be observed with the naked eye, including 
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some forms of mirabilite that can be rather large in size (Fig. 6d) and some intermediate states 
between mirabilite and anhydrous sodium sulfate that can exhibit grainy or layered structures 
(Figs. 6e and f). 
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(a) Mirabilite in soils (b) Anhydrous sodium sulfate crystals 


(c) Powdered sodium sulfate (d) Mirabilite in solution 


(e) Granulated sodium sulfate (f) Layered sodium sulfate crystals 


Fig. 6 Crystal morphology of sodium sulfate 


3.2 Effect of chloride ion content on salt expansion 
3.2.1 Supersaturation ratio in Na?28O4-NaCI binary solutions 


Sodium chloride crystals cannot be formed, as they show consistently high solubility in various 
sodium chloride solutions. Nevertheless, the sodium chloride content has a significant effect on 
the solubility of sodium sulfate. The solubilities of sodium sulfate solutions under different 
sodium chloride contents and different temperatures are given in Table 1. 

The saturation concentration of a sodium sulfate solution gradually decreases with increasing 
concentration of sodium chloride. When the sodium chloride concentration exceeds 6 mol/L, the 
saturation concentration of the sodium sulfate solution is very low. Consequently, 6 mol/L of 
sodium chloride was selected as the maximum concentration for the sodium chloride solution. 
Then the water activity and mean ion activity of the sodium sulfate solution in a Na?S5O4-NaCI 
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binary solution were calculated from Equation 1, and the related parameters are listed in 
Appendix B. Given the cooling rates are same, the calculated supersaturation ratios of the sodium 
sulfate solution are shown in Figure 7. 


Table 1 Saturation concentration (mol/L) of Na»SO4 under different sodium chloride contents and different 
temperatures 


5?C 10°C 15°C 21.5°C 21°C 35°C 


NaCl NaSO, NaCl NaSO, NaCl Na3SO, NaCl NaSO, NaCl NaSO, NaCl NaSO, 


0.20 0.42 0.00 0.64 0.00 0.96 0.00 1.50 0.00 2.18 0.00 3.38 
0.86 0.32 0.79 0.44 0.77 0.67 1.62 1.09 0.48 2.02 0.38 3.08 


1.52 0.25 1.76 0.33 1.72 0.50 3.13 0.97 0.95 1.91 2.43 1.88 
2.65 0.20 2.80 0.28 2.80 0.42 3.66 0.96 1.42 1.83 3.37 1.39 
3.53 0.19 3.91 0.28 3.91 0.42 4.66 1.06 2.89 1.75 5.72 0.58 
4.52 0.21 4.96 0.30 5.04 0.44 4.75 1.02 3.39 1.49 6.39 0.00 
5.05 0.22 5.40 0.31 5.40 0.46 5.70 0.72 3.52 1.42 - - 
5.87 0.22 5.80 0.32 5.78 0.48 6.04 0.31 3.72 1.36 - - 
6.08 0.20 6.08 0.23 6.05 0.25 6.35 0.00 5.79 0.65 - - 
6.50 0.00 6.45 0.00 6.35 0.00 - - 
Note: “—” indicates no value. 


From Figures 7 a and b, the supersaturation ratios of the sodium sulfate solution decrease when 
the concentration of the sodium chloride solution is less than 2.0 mol/L under 27°C and 21.5°C 
conditions, and inflection points are observed for a 5.0 mol/L concentration of sodium chloride 
solution. Above this concentration, the supersaturation ratios of the sodium sulfate solution 
markedly increase as the concentration of the sodium chloride solution increases. However, when 
the concentration of the sodium sulfate solution exceeds 2.0 mol/L, a critical point occurs at the 
3.0 mol/L sodium chloride solution. In other words, when the concentration of the sodium 
chloride solution is lower or higher than 3.0 mol/L, the supersaturation ratios of the sodium 
sulfate solution decrease or increase correspondingly (Figs. 7a and b). 

From Figures 7c-e, variations in supersaturation ratios of the sodium sulfate solution are 
similar, i.e., they increase as the concentration of sodium chloride solution is lower than 3.0 mol/L 
and then decrease until the concentration reaches up to 5.5 mol/L. Furthermore, there is an abrupt 
increase in the supersaturation ratio as the concentration is higher than 5.5 mol/L. On the other 
hand, the supersaturation ratios increase as the concentration of sodium sulfate solution is less 
than 1.5 mol/L. As shown in Figures 7a—e, the maximum supersaturation ratio at 5°C is about 40 
times higher than that of at 27?C. Thus, we speculated that the supersaturation ratios will 
constantly increase at a given concentration with decreasing temperature (the concentration 
corresponding to binary system solutions). 


3.2.2 Effect of sodium chloride content on salt expansion 


The salt expansion ratio of sulfate saline soils increases concurrently with the increase in sodium 
chloride content (Fig. 8). Salt expansion at temperatures above 20°C only occurs in the sample of 
sulfate saline soil containing 3.8% sodium chloride. This is in accordance with the rule for 
supersaturation ratio of sodium sulfate solution, presented in the Na»SO4-NaCI binary solution. 
The supersaturation ratio of the sodium sulfate solution increases as the chloride concentration 
increases. When this ratio exceeds the initial supersaturation ratio, mirabilite precipitates and 
causes salt expansion. However, in the sulfate saline soil containing 5.8% of sodium chloride, the 
sodium sulfate solution is saturated. In this case, the sodium sulfate exists in the form of a powder 
in soils above 20°C (in Fig. 6c) and soil displacement remains constant. With declining 
temperature, more powdered sodium sulfate can dissolve in water to form mirabilite. At this point, 
due to the greater precipitation of mirabilite, the salt expansion ratio increases sharply and 
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Fig. 7  Supersaturation ratios of NazSO4 in Na?SO4-NaCI binary solution under different temperatures 
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Fig.8 Salt expansion of sulfate saline soil under different sodium chloride contents 
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4 Discussion 


The behavior of salt solution in soils is usually affected by the temperature, pressure, and 
concentration. The soil pore pressure will change with salt crystallization, but this has little effect 
on the equilibria between liquid and solids. We therefore do not need to consider the effect of 
pressure on crystallization. The results present above suggest that the supersaturation ratio of a 
sodium sulfate solution is strongly associated with crystal morphology, cooling rate, and chloride 
ion concentration. It appears that a small supersaturation ratio generally leads to stable crystal 
shape, and that large unsteady-state zones are associated with the formation of unstable crystals, 
such as sodium sulfate heptahydrate. The cooling rate not only affects the supersaturation ratio of 
the sodium sulfate solution, but also influences the extent of the non-steady-state zone. An 
increase in cooling rate likely results in the formation of unstable crystal shapes, and may reduce 
salt expansion. 

The temperature of the precipitated salt crystals in soils is generally lower than that in pure 
solution (Wan and Lai, 2013), and the adsorptive action of free energy on soil particle surfaces 
leads to increased supersaturation ratios for salt solutions in soils. This, in turn, results in an 
enlarged non-steady-state zone for the sodium sulfate solution, showing that salt crystallization is 
a more complex process in soils than in pure solutions. Crystallization theory states that salt 
crystals are first generated in the larger pores of porous materials and that crystallization becomes 
more difficult with the decreasing pore size (Steiger, 2005). During salt expansion, crystals 
continue to fill the pores and decrease the pore volume, so the salt expansion ratio increases 
slowly with time (Fig. 4). When the temperature is below the freezing point of saline soils, water 
begins to freeze and growth of ice crystals decreases pore size. On the other hand, water freezing 
increases the concentration of the remaining salt solution. Thus, these two aspects result in a high 
initial supersaturation ratio (Fig. 1) and thus a weak salt expansion. 

The condensation depression of the solution decreases as cooling rate decreases, which can 
promote the precipitation of salt crystals. The mass of crystals increases more quickly in the 
slowly cooled soils than in the rapidly cooled soils. In addition, because the volume of one mole 
of sodium sulfate decahydrate is 3.1 times greater than that of an anhydrous sodium sulfate (Chen 
et al., 2006), the transformation of anhydrous sodium sulfate into mirabilite enhances the salt 
expansion process. The greatest variation of salt expansion occurs as the precipitation of salt 
crystals begins, and the concentration of salt solution is close to the saturation concentration for 
crystallization. We can therefore estimate the saline soil deformation in saturated soil during 
cooling based on the actual solubility curve, simply by subtracting the volume of water loss from 
the increase in the volume of salt crystals. It remains very difficult to calculate the amount of salt 
heaving in unsaturated soil. 

The presence of sodium chloride influences both the water and ionic activities of the sodium 
sulfate solution. Indeed, with increasing concentration of the sodium chloride solution, the 
maximum temperature for mirabilite to precipitate first decreases from 32.4?C to 17.9?C, and then 
increases for supersaturated concentrations. In solutions with high salt content (concentrations of 
the sodium sulfate solution exceed than 2.0 mol/L), the increased chloride ion content can 
effectively reduce the supersaturation ratio of the sodium sulfate solution at temperatures above 
17.9?C, thereby effectively inhibiting salt expansion. However, below 17.9?C, the mirabilite 
exists in a stable state and sulfate saline soils readily undergo salt expansion. 

In seasonally frozen regions, the temperature is sometimes below zero. In these settings, an 
increase in chloride ion content may lead to the crystallization of sodium chloride dehydrate and 
cause salt expansion. Consequently, it is very important to use chloride ion content in a 
reasonable manner, considering the temperature range of the intended environment, to control the 
strength of salt expansion and frost heaving in sulfate saline soils. In addition, increasing the dry 
density of subgrade fill and decreasing the median pore size is an effective way to weaken salt 
expansion. 
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5 Conclusions 


Based on the cooling experiments and theoretical analyses of sulfate saline soils, we found that 
the initial supersaturation ratio of the brine in sulfate soils increases linearly as temperature 
decreases, and this increase is more rapid below the freezing temperature of the soil. The non 
steady-state zone of sulfate saline soils is determined by the cooling rate; the higher the cooling 
rate, the larger the initial supersaturation ratio of the sodium sulfate solution. The supersaturation 
ratio of the sodium sulfate solution is closely related to the crystal morphology of the 
precipitating mineral phases. The initial supersaturation ratio decreases as the stability of the 
crystal shape increases. When salt crystals precipitate as mirabilite, the salt expansion ratio is at 
its maximum and the effect of salt expansion is significant. Chloride ion concentration has a 
distinct influence on the crystallization supersaturation ratio of the sodium sulfate solution. Below 
17.9°C, an increase in chloride ion content induces soil heaving in sulfate saline soils. 
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Appendix A 
The detailed definition of water activity in solution can be found in Prausnitz et al. (1999). It is 
related to the osmotic coefficient o (Steiger, 2005), the value can be calculated by Equation A1: 
Lo M» Al 
Ina, = "91995 24^: ad 
where M,, is the molar mass of water, and m; is the molar mass fraction of i™ ion. In the Pitzer 
approach, the osmotic coefficient for a single electrolyte is given by: 


3/2 
p-1=|z,z |f «2m "Bt, «one 992 os, (A2) 
v v 


where f? =A 


, and Bj, = Bry + Bux expla D) + yy exp(a I): z«z-is the charge number 
1+bI 


of electrolyte ions; m is the molar mass fraction of the solute; v is the total number of ions 
released upon complete dissociation of the salt; and v+v- is the number of positive and negative 


ions; b is a constant (b=1.2 kg!?mol!?); and a, and az are regression coefficients with different 


values for different ionic types. The interaction parameters 8? , 40 , B° and C$, are 


specific to a salt MX, and are a function of the temperature. The relation between the interaction 


parameters 8% , B) and Cf,,and temperature is expressed by the following equation: 


bos : 
P(T)=a,+@,0/T-1/T,) - a,ln(T/T,) - (T -T,) - (T^ — T; ) + a, ln(T —225), (A3) 
where a\—d¢ refer to the regression parameters (Steiger et al., 2008). Their values vary in different 
solutes. Tr=298.15 K. 
The mean ion activity is defined as: 
a, — uy. (A4) 
where m+ is the mean mass molarity; and y+ is the mean ionic activity coefficient. 
For 2-1 salt types, such as Na2SOa, m+ can be approximated as (Prausnitz et al., 1999): 


13 
m, =4 ` Myx > 


(A5) 
where MX is the corresponding salt. 
The mean ionic activity coefficient can be calculated using the following equations (Marliacy, 


2000): 


3/2 
In y, =| f + 2m Blag, ease Dt cy, (A6) 
where, 
1/2 2 " 
"=A,(——, t -In-bl ^), 
Pour ig OPEM 


D 
212/824 2s al (1 a1? -0.5a71)xexp(— -a,P^) |+ Tx 


[1- q^ i —-0.5a;1) |xexp(-a,1'”), 


3 
Cix ES 2 Chy . 


The Pitzer model is used to calculate the supersaturation ratio at positive temperatures (Pitzer, 
1975; Steiger, et al., 2008), while the FREZCHEM or SMW (Structural Modeling Workflow) 
model is used to calculate the supersaturation ratio below the freezing point of pure water 
(Spencer et al., 1990; Marion and Grant, 1994; Marion and Farren, 1999). 


A, is the Debye-Hückel parameter for the osmotic coefficient, which represents a nonlinear 


relationship with the temperature T (Clegg and Whitfield, 1991): 
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205.381 (A7) 


A, = 0.04944T —1.963In T — 10.792. 


Iis the ionic strength defined as: 
-lymg. (A8) 
2*1 
where z; is the charge number of electrolyte ions i. 


Appendix B 


For a Na2SOu-NaCl binary solution system, the osmotic coefficient g can be calculated using the 
formula defined by Marliacy et al. (2000). The mean molarity (m+) and the mean ion activity (y+) 
coefficients are respectively defined as: 
m, =(m} m>)"; r=)» (B1) 

where m+ and m refer to the molar mass of positive and negative ions, respectively; and y+ and y_ 
refer to the ion activity coefficient of positive and negative ions, respectively. The remaining 
parameters are defined above. 

The ion activity coefficient of sodium and sulfate ions can be calculated by Equations B2 and 
B3, Ea 


Fps m2 X F42m (Bua + ZC ya) + Ugo. B, s ZC ae) 


cr NaCl 

, (B2) 

x LX " x x X 4X 

TUE Mso} Pia cr SOP Fli na ” Mer Cw) TM Mso} Cwsso, 
z 

= x X x x 

In ^ SO; < sor F+ 2m, a* (Byaso, T ZC, so, p 2m. Di „SOF * Mop xm. Todd (B3) 
, 


+Z soz | Mya xm. X Crac) m, . X Mo} X Cy, so, 


M 


where Bmx is the second virial coefficient; Cyx is the third virial coefficient; ^ — and 
CT ,SO; 


sor 


Pc sox are the fitting parameters for the binary solution of sodium sulfate and sodium chloride, 
Na*,Cl ,SO;- 


respectively; Bmx and Cyx are calculated by Equations B4 and BS, respectively (Marliacy et al., 
2000; Steiger et al., 2008): 


By = BO. + Bia e (o D) Bi e (o1), (B4) 
NE^ (B5) 
MX y2"* 
2|z,2;| 


The parameter Qo sow Can be calculated as: 
Oz 


ó = E Eg 
Cl.SOP — Oor sor T. Oor sor *Ix E ? (B 6) 
where "0, "s is the term accounting for asymmetrical mixing (Pitzer, 1975); and ^9 is 
SUA as 
the derivative of ^9 


CI-,S042- * 


The function F in Equations B2 and B3 is defined as: 


2 


--A PET are += nu bI ^) em, _x Bug tim, x Dos. X By, so, + Mo X Mop. X D so . (B7) 
with p. =A. D o (a1) + BO g (a - and o sor = zd E where g' is the derivative of g. 


The "ew g in Equation B4 is defined as: 
g-2|1-0*3e7 x^. (B8) 


I". The remaining parameters in Equations B2 and B3 are defined as above. 


with x =a)!” or a2 


